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Abstract 

in 

Considering the hypothesis of mixing of three active neutrinos with, at least, one 



sterile neutrino, we report on a simple 4x4 texture whose 3x3 part arises from the 
popular bimaximal texture for three active neutrinos v e , v^, v T , where C12 = l/v2 = S12, 



O 

c 23 = l/v2 = s 2 3 and Si 3 = 0. Such a 3 x 3 bimaximal texture is perturbed through a 



rotation in the 14 plane, where v± is the extra neutrino mass state induced by the sterile 
neutrino v s which becomes responsible for the LSND effect. Then, with m\ ~ m\ we 
predict that sin 2 20 atm = ~(1 + cf 4 ) ~ 0.99 and sm 2 26- LSN D = \s\ A ~ 4.5 x 10" 4 , and 
in addition Am^ tm = Am| 2 and AmL SND = (Am^l, where cf 4 = sin 2 26 l so i ~ 0.97 and 
Am^ = Am 2 ol ~ 10~ 7 eV 2 if e.g. the LOW solar solution is applied. In this four-neutrino 
texture with m\ ~ m\ the sum rule sin 2 26* S oi + \ sin 2 26 l chooz + sin 2 26 i lsnd = 1 holds in 
the two-flavor approximation (for each of three cases), leaving room for the LSND effect, 
depending on the magnitude of Chooz effect that, not observed so far, leads (at present) 
to the upper bound sin 26*lsnd ~ 1.3 x 10 3 and the lower bound sin 2 2# so i ~ 0.95. At 
the end a four-neutrino seesaw mechanism is sketched. 
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The present status of experimental data for atmospheric v^s as well as solar z/ e 's 
favours oscillations between three conventional neutrinos v e , , v T only [1]. However, the 
problem of the third neutrino mass-square difference, related to the possible LSND effect 
for accelerator z/^'s, is still actual [2], stimulating a further discussion about mixing of these 
three active neutrinos with, at least, one hypothetical sterile neutrino v s (although such a 
sterile neutrino is not necessarily able to explain the LSND effect [3]). As a contribution 
to this discussion, we report in this note on a simple 4x4 texture for three active and one 
sterile neutrinos, v e , , v T and is s , whose 3x3 part arises from the popular bimaximal 
texture [4] working grosso modo in a satisfactory way for solar z/ e 's and atmospheric z/^'s 
if the LSND effect is ignored. Such a 3 x 3 bimaximal texture is perturbed [5] by the 
sterile neutrino v s inducing one extra neutrino mass state v± and so, becoming responsible 
for the possible LSND effect. In fact, with the use of our 4x4 texture we predict that 
sin 2 26 l LSND = \ s u and Am 2 SND = (Amf-J, while sin 2 2# so i = c 2 4 and Am 2 ol = Am^ as 
well as sin 2 26' atm = |(1 + c 2 4 ) and Am 2 tm = Am| 2 , if m\ ~ m\ (and both are different 
enough from m\ and m 4 ). Here, c 2 4 ~ 0.97 and Am^ ~ lCT 7 eV 2 if e.g. the LOW solar 
solution [1,6,7] in its recent estimation [7] is accepted; then we predict sin 2 2# atm ~ 0.99 
and sin 2 2# LSN d ~ 4.5 x 10" 4 . 

In the popular 3x3 bimaximal texture the mixing matrix has the form [4] 

/ l/y/2 l/y/2 \ 

C/( 3 ) = -1/2 1/2 l/y/2 . (1) 

V 1/2 -1/2 LA/2/ 
Such a form corresponds to c i2 = 1/a/2 = Si 2 , c 23 = 1/V2 = s 2 3 and s 13 = in the 
notation used for a generic Cabibbo-Kobayashi-Maskawa-type matrix [8] (if the LSND 
effect is ignored, the upper bound | s 13 1 ~ 0.1 follows from the negative result of Chooz 
reactor experiment [9]). Going out from the form (1), we propose in the 4x4 texture the 
following mixing matrix: 



/ 0\ / c u s u \ 

10 

1 

V 1 / V -su cu / 



/ c l4 /V2 l/y/2 s u /V2\ 

-Ci 4 /2 1/2 1/V2 -su/2 

c M /2 -1/2 l/y/2 Su/2 
\ su cu J 



1 



with C14 = cos 9 14 and s i4 = sin 6*14 (note that in Eq. (2) only Si 2 , s 23 and s i4 of all 
with i,j — 1, 2, 3, 4 , 2 < j are nonzero). The unitary transformation describing the mixing 
of four neutrinos v a = v e , , z/ r , z/ s is inverse to the form 

"a = Yl U »i V i » ( 3 ) 

where i/j = z/i , z/ 2 , z/ 3 , z/ 4 denote four massive neutrino states carrying the masses m 8 = 
mi , m 2 , m 3 , m 4 . Here, [/ = (£/ Q i) , a = e , /i , r , s and « = 1 , 2, 3, 4. Of course, W = 
U^ 1 and also U* = U, so that a tiny CP violation is ignored. 

In the representation, where the mass matrix of three charged leptons e" , ji~ , r~ is 
diagonal, the 4x4 neutrino mixing matrix U is at the same time the diagonalizing matrix 
for the 4x4 neutrino mass matrix M = (M a p): 

U'MU = diag(mi , m 2 , m 3 , m 4 ) , (4) 

where by definition m\ < m\ < m\ and either m\ < m\ or m\ < m\. Then, due to the 
formula M a/3 = J2i UaiiriiU^ we obtain 



M, 



(if! 



M, 



fifi 



i (c? 4 mi + s 2 u rri4 + m 2 ) , 
-M er = ~^/= (ci 4 mi + Si 4 m 4 - m 2 ) , 
cf 4 mi + Si 4 m 4 + m 2 ) + m 3 



M, 



(IT 



M P 



M^ s = 
M xx = 



i (c^rri! + Si 4 m 4 + m 2 ) - m 3 



-^|Ci4Si 4 (mi -m 4 ) , 

-M TS = i C14 S14 (mi - m 4 ) = --^=M es , 
Si 4 mi + c 2 14 m 4 . 



(5) 



Of course, M M and also M* = M. From Eqs. (5) we find that 



M ee - M e ^ + M ss 
mi,4 or m 4 ,i = - ± 



( M ee - M efl V2 - M s . 



m 2 



= M ee + M eix V2 , m 3 = M^ + M^ 



+ 2M*. , 



(6) 



if 777.4 < m i or 777 1 < 7774, respectively, and 



( 2c S ? = 8M e 2 * (J) 

14 1 (M ee - M e ,V2 - M ss f + 8MI ' 1 ' 

Obviously, mi + m 2 + m 3 + m 4 = M ee + M w + M TT + M ss as M ee = M w — M MT and 

M w = M rr . 

Due to the mixing of four neutrino fields described in Eq. (3), neutrino states mix 
according to the form 

k> = £^h>- (8) 

i 

This implies the following familiar formulae for probabilities of neutrino oscillations v a — > 
Ufs on the energy shell: 

P (y a -+v p ) = \ {v p \e iPL \v a )\ 2 = 8 Pa - 4 £ U^U^UojU^ sin 2 ^ , (9) 

valid if the quartic product 1/^1/^1/^1/^ is real, what is certainly true when a possible 
CP violation can be ignored (then U* { = U ai ). Here, 



Am 2 -L 
~E~ 

with Am 2 j, L and measured in eV 2 , km and GeV, respectively (L and E denote the 



Xji = 1.27 — , Am 2 . = m 2 - m 2 (10) 



experimental baseline and neutrino energy, while Pi = \J E 2 — rri 2 ~ E — m 2 /2E are 
eigenstates of the neutrino momentum P). 

With the use of oscillation formulae (9), the proposal (2) for the 4x4 neutrino mixing 
matrix leads to the probabilities 



P(v e — > v e ) ~ 1 - c 2 4 sin 2 x 2 i - (l + c 2 4 ) s 2 4 sin 2 x 44 , 

P{y^i -> ^v) = -> f T ) — 1 - ^c 2 4 sin 2 x 2 i - ^(1 + c 2 4 ) (^sin 2 x 32 + ^s 2 4 sin 2 x 41 ^ 



^s 2 4 sin 2 x 43 , 
- — > i/„l ~ 

2 



P{.Vfi —> v e ) = P(v T -> z/ e ) ~ i (c 2 4 sin 2 x 2 i + sin 2 x 4 i) , 
f(^ -> f T ) - -^c 2 4 sin 2 a;2i + ^(1 + c 2 4 ) ^sin 2 x 32 - ^s 2 4 sin 2 a; 4 i^ 



+ is 2 4 sin 2 x 43 (11) 
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in the approximation, where m 2 ~ m 2 , (and both are different enough from m| and m|). 
The probabilities involving the sterile neutrino v s read: 

-> = -> v a ) = (C14S14) 2 sin 2 x 4 i , 

P(v e ^v s ) = 2(ci 4 si 4 ) 2 sin 2 a;4i , 

P(v s -> i/ s ) = 1 - (2ci4Si 4 ) 2 sin 2 x 4 i . (12) 
If Am 2 .! < lAmlil {i.e., x 21 < |ar 41 |) and [1,6,7] 

A?t4 = Am s 2 ol ~ (1(T 5 or 10~ 7 or 10~ 10 ) eV 2 , (13) 
then, under the conditions of solar experiments the first Eq. (11) gives 

P(v e -> z/ e ) sol ~ 1 - c 44 sin 2 (aj2i) so i - ^(1 + c 2 4 )s 2 4 (14) 

with the estimate 

c 2 4 = sin 2 2d sol ~ (0.66 or 0.97 or 0.80) , ^(1 + c 2 4 )s 2 4 ~ (0.28 or 0.030 or 0.18) , (15) 

if the constant term in Eq. (14) can be considered as a small perturbation of the usual 
two-flavor formula. In Eqs. (13) and (15) the recent estimation [7] for the LMA or LOW 
or VAC solar solution, respectively, is used. Note that 

P(y e -> z/ e ) so i ~ 1 - P(y e -> z/ M ) sol - P(y e -> tv)soi - (ci 4 Si 4 ) 2 (16) 

with (c 14 s 14 ) 2 ~ (0.22 or 0.029 or 0.16). 

If Am 2 x < Am 2 2 < (Am^l , |Am| 3 ] (i.e., x 2 i < ^32 < |^4i| , |#43|) and [1] 

Am 2 2 = Am 2 tm ~ 3 x lO" 3 eV 2 , (17) 
then for atmospheric experiments the second Eq. (11) leads to 

P{ypL -> ^)atm ~ 1 - ^(1 + c 44 ) sin 2 (a;3 2 ) a t m - ^(3 + c 2 4 )s 2 4 (18) 
with the prediction 
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sin 2 26 atm = 1(1 + c 2 4 ) ~ (0.83 or 0.99 or 0.90) i(3 + c 2 4 )s 2 4 ~ (0.16 or 0.015 or 0.095) 
2 8 

(19) 

following from the value (15) of c 2 4 , if again the constant term in Eq. (18) can be 
considered as a small perturbation. Notice that 

P{y? -»• ^)atm ^ 1 - PiVp -> ^r)atm ~ ^(1 + C 2 4 )s? 4 (20) 

with (1 + c 2 4 )sf 4 /4 ~ (0.14 or 0.015 or 0.09). 

Eventually, if Am^ <C |Am 41 | {i.e., X21 <C |#4i|) and [2] 

|Am 2 4 | = Am 2 SND > 0.1 , e.g. ~ 1 eV 2 , (21) 
then for the LSND accelerator experiment the third Eq. (11) implies 

P(z/ M -> Z/ e ) LS ND ^ ~sf 4 sin 2 (X41) LSND (22) 

with the prediction 

sin 2 2# LS nd = -Su ~ (0.058 or 4.5 x 10~ 4 or 0.020) (23) 
inferred from the value (15) of c\ A . Note that 

1 / s \ ^ 

P(z/ M -> z/ 6 )lsnd ^ - (^J Piy -> za*)lsnd (24) 

with i(si 4 /ci 4 ) 2 ~ (0.13 or 4.8 x 10~ 4 or 0.031). 

It may be worthwhile to remark that if there were c 2 4 = sin 2 2# so i ~ 0.9 (as in the case 
of older estimations for the LOW solar solution [1]), our predictions would be sin 2 2# atm = 
|(1 + c 2 4 ) ~ 0.95 and sin 2 2# LSN d = |s 4 4 ~5x lO" 3 . 

Concluding, we can say that Eqs. (14), (18) and (22) are not inconsistent with so- 
lar, atmospheric and LSND experiments. All three depend on one common correlating 
parameter c 2 4 , implying c 2 4 = sin 2 1Q so \ ~ (0.66or0.97or0.80), sin 2 2# atm = \{\ + c 2 4 ) ~ 
(0.83 or 0.99 or 0.90) and sin 2 2# LSN d = §4i ~ (0.058 or 4.5 x 10~ 4 or 0.020). They depend 
also on three different mass-square scales Am^ = Am s 2 ol ~ (10~ 5 or 10 -7 or 10~ 10 ) eV 2 , 
Amj 2 = Am 2 tm ~ 3 x 10~ 3 eV 2 and (Am 2 ^ = AmJ SND ~ 1 eV 2 . Here, the LMA or LOW 
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or VAC solar solution [7] is accepted (perturbatively). Note that in Eqs. (14) and (18) 
there are constant terms which modify moderately the usual two-flavor formulae. Any 
LSND-type accelerator project, in contrast to the solar and atmospheric experiments, 
investigates a small — > v e oscillation effect caused possibly by the sterile neutrino v s . 
Thus, this effect (if it exists) plays the role of a small "sterile" perturbation of the basic 
bimaximal texture for three active neutrinos v e , , v T . Of course, if su were zero, such 
an LSND effect would not exist and both solar u e — > u e and atmospheric — > z/ M oscil- 
lations would be maximal. So, from the standpoint of our texture (2), the estimated not 
full maximality of solar v e — > u e oscillations may be considered as an argument for the 
existence of the LSND effect. 

The final results (14), (18) and (22) follow from the first three oscillation formulae 
(11), if either 

m 2 -C m 2 ~ m 2 . ~ m 2 (25) 

with 

m\ ~ 1 eV 2 , m\ < 1 eV 2 , Am 2 21 ~ (1(T 5 or 1(T 7 or 1(T 10 ) eV 2 , Am 2 2 ~ 3 x 1(T 3 eV 2 

(26) 

or 

ml ~ ml < mj < m 2 (27) 

with 

ml < 1 eV 2 , mj ~ 1 eV 2 , Am 2 21 ~ (10^ 5 or 1(T 7 or 10^ 10 ) eV 2 , Am 2 2 ~3x 10 -3 eV 2 . 

(28) 

In both cases Am| x Am 2 2 (Am^j ~ 1 eV 2 . The first case of m\ <C m 2 ~ 1 eV 2 , 
where the neutrino mass state z/ 4 induced by the sterile neutrino v s gets a vanishing mass, 
seems to be more natural than the second case of m\ <^ ml ~ 1 eV 2 , where such a 
state gains a considerable amount of mass ~ 1 eV "for nothing". This is so, unless one 
believes in the liberal maxim "whatever is not forbidden is allowed". Note that in the first 
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case the neutrino mass states v\ , v 2 , v$ get their considerable masses ~ 1 eV through 
spontaneously breaking the electroweak SU(2)^xU(l) symmetry which, if it were not 
broken, would forbid these masses. 

Finally, for the Chooz reactor experiment [9], where it happens that (^ji)chooz — 
(^ji)atm, the first Eq. (11) predicts 

P(v e -> V e )chooz ^ Pipe -> ^e)atm ~ 1 - ^(1 + c\ A )s\ A (29) 

with |(1 + c\ A )s\ A ~ (0.28 or 0.030 or 0.18) from Eq. (15). Here, sin 2 (:r 41 )chooz = \ since 

|(z4l)chooz| - |(^4l)atm| > (^atm ~ 1 with |Am 41 | > Affl^. 

In terms of the usual two-flavor formula, the negative result of Chooz experiment ex- 
cludes the disappearance process of reactor P e 's for moving (l+c 2 4 )s 2 4 = sin 2 2#cho OZ ~0.1, 
when the range of moving |Am 41 | = AmQ hooz ~ 0.1 eV 2 is considered ( then AmQ hooz ^> 
Am 2 tm ~ 3 x 10~ 3 eV 2 , implying sin 2 xchooz = §)■ Thus, the nonobservation of Chooz ef- 
fect for reactor z/ e 's in the above parameter ranges leads to (1 + c 2 4 )sf 4 ~ 0.1 and hence, to 
the upper bound sin 2 2#lsnd = \ s \a ~ 1-3 x 10 -3 , when AmL SND = |Am 41 | = Am^ ~ 
0.1 eV 2 . It means that sin 2#lsnd, constrained by Chooz (in our four-neutrino texture), 
lies outside the parameter region suggested at 90% CL by the existing LSND data [2], 
if the KARMEN2 results [2] excluding a large part of this region are taken into account 
(in fact, in this corrected region sin 2 26*lsnd ~ 2 x 10~ 3 ). But, at 99% CL, this may be 
not true, allowing for Am^sND ~ 1 e ^ 2 ( as i then, in the existing LSND parameter region 
sin 2 2#LSND ~ 8 x 10~ 4 ). At any rate, among three solar neutrino solutions considered 
here [7], only the LOW solution is consistent with the Chooz bound [cf. Eq. (23)]. Also 
the value sin 2 26*lsnd ~ 5 x 10~ 3 , mentioned before as corresponding to older estimations 
for LOW solar solution [1]: sin 2 26 so \ ~ 0.9, is eliminated by this bound. 

However, the Chooz-allowed, LOW-induced value sin 2 2#lsnd ~ 4.5 x 10~ 4 corre- 
sponding to sin 2 26' so i ~ 0.97 is situated (in contrast to sin 2 xlsnd ~ 5 x 10~ 3 ) outside the 
parameter region implied by the existing LSND data [2] (even at 99% CL). 

Of course, the existence of Chooz bound for the LSND effect and of the relation of the 
latter to solar neutrino solutions is caused by the correlations between different neutrino- 
oscillation modes connected through the parameter s 2 4 appearing in our four-neutrino 
texture [cf. Eqs. (11)]. In fact, the identities |s 4 4 = |(1 — c 2 4 ) 2 , (1 + c 2 4 )s 2 4 = 1 — c 4 4 
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and c 2 4 + |(1 + c\ A )s\ A + \s\ A = 1 can be translated into the correlations 

sin 2 2# LSND = - (1 - sin 2 26 sol ) 2 , sin 2 2# Chooz = 1 - sin 4 26 sol (30) 
and the sum rule 

sin 2 29 sol + ^ sin 2 29 Chooz + sin 2 2# LSN d = 1 (31) 

for three neutrino-oscillation amplitudes (each in the reasonable two-flavor approxima- 
tion). The sum rule (31) can be derived also from the probability summation relation 
J2p P(y e ~^ v p) — 1 (with j3 — e , /i , t , s) considered under the assumption m 2 ~ m\ for 
solar z/ e 's (when |(x 4 i) so i| > (x 2 i) so i - 7r/2). 

We can see that, when accepting the present Chooz results, we stand with our four- 
neutrino texture before the alternative: either there is no LSND effect at all (then 
sin 2 26 l LSND = \ s \a = and we are left with the three-neutrino bimaximal texture [4,8]), 
or this effect exists all right, but at a point in the parameter space, where the oscillation 
amplitude sin 2 2#lsnd is shifted (versus the existing LSND data) towards a smaller value 
sin 2 26 l LSND = \s\ A ~ 1.3 x 10" 3 (though >0). Note that, if sin 2 2# LS nd = \s\ A was at 
the Chooz bound value 1.3 x 10~ 3 , then sin 2 26 so \ = c\ A would be at 0.95. If, rather, 
sin 2 2# LSND = \s\ A was at the value 8 x 10~ 4 equal to its existing LSND lower limit at 99% 
CL, then sin 2 26 so \ = c 2 4 would be at 0.96. Of course, we should keep in mind the fact 
that the present estimates for sin 2 29 so \ (and even more for sin 2 26 l LSND) are preliminary 
and, in fact, very fragile. 

From the neutrinoless double (3 decay, not observed so far, the experimental bound 
M ee = | U^rrii] ~ [0.4 (0.2) - 1.0 (0.6)] eV follows [10] (here, C/ 2 appears even if U* ei ^ 
U e i). On the other hand, for e.g. the LOW solar solution the first Eq. (5) gives 

— , , 1. 

M ee = \M ee \ ~ -|0.97mi + 0.03m 4 + m 2 | , (32) 

what in the case of Eq. (25) with mi~±l eV and m 2 ~l eV or Eq. (27) with |m 4 |~l eV 
leads to the estimation M ee ~ (0.99, 0.015) eV or M ee ~ 0,015 eV, respectively (putting 
M ee = \M ee \ in Eq. (32) we ignore a possible violation of CP: we get [/* = U ei , since 
M ee = E i \U ei \ 2 m i ). 
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In the case of neutrinoless double (3 decay we assume, of course, that in our texture 
four flavor neutrinos v a — v e , , v T , v s are Majorana fermions v a = v aL + (zvl) C , where 
V&R = {yah) c = (O-R an d v & — u a- Then, the neutrino mass term in the Lagrangian 
density, 

- 4K? = \H^ M aP (33) 

Z a/3 

with M = {M a p) given in Eq. (5), is a Majorana lefthanded mass term, as v a are built 
up of v aL [and [y aL ) c \. 

It may happen, however, that in this case, beside four flavor neutrinos v a correspond- 
ing to four light mass neutrino states Vi = v\ , u 2 , 1/3, z/ 4 , there exist four other flavour 
neutrinos v' a — v' e , 1/ , v' T , i/, being also Majorana fermions u' a = u' aR + (v' aR ) c , where 
v' aL = {v' aR ) c = l and v' a = iA, but connected this time with four heavy mass neutrino 
states v[ — v[ , v' 2 , v' z , 1/4. Thus, the latter may be practically decoupled from the former 
light mass neutrino states. In other words, a jour-neutrino seesaw mechanism may work, 
"explaining" the lightness of z/, versus u[: 



M aj3 = M% sht) ~ M$ - (M^M^-'M^) ~- (M^^-'M^l) , 

ee ~ Mif , (34) 

where M^ L ^ <C M^ D ^ <C M*^ in the perturbative sense. These 4x4 matrices, viz. 
Majorana lefthanded, Dirac and Majorana righthanded, respectively, are assumed to be 
real and symmetric for simplicity. Here, the overall neutrino mass term in the Lagrangian 
density has the form 

_r = - Vfz/ d') I M °A W up \ £ (light) - r( heav y) 

2 a/3 V M «/3 ) \ P / 

Note that the combinations = v a L + 1/ R are formed in the sum of two terms with 
M^jP in Eq. (35) and play there the role of four Dirac flavor neutrinos. 

For three conventional a = e , fj, , r the Majorana neutrinos u a = u a i + [yahf = 
vL D l + (^JlY and 1/ = 1/^ + (u' aR ) c = u ( a D R ] + (u { a D R ] ) c were denoted also by i/(°> and z/W, 
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and called active and conventional-sterile neutrinos, while the Dirac neutrinos were 
written simply as v a (cf. e.g. Ref. [5]). Then, v' aR = were written as v aR . 

The sterile neutrino v s — the extra flavor neutrino considered in this paper beside 
the active neutrinos v e , , v T — differs in an obvious way from any of four possible 
conventional-sterile neutrinos v ' e , 1/ , v' T and v' s needed for the four-neutrino seesaw mech- 
anism to be realized (and for four Dirac neutrinos , , and to be defined) . 
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Appendix: Allowing for nonzero Si 3 

If we maintain two maximal mixings C\ 2 = 1/a/2 = s 12 and C23 = 1/V2 = s 2 3, but 
decide to allow Si 3 to be nonzero (though small, as may be expected), then Eqs. (1) and 
(2) take the forms 



C/(3) 



/ c 13 /V2 

-(l + s 13 )/2 
V (l-Si 3 )/2 



C^/v 7 ^ S13 \ 

(l-s 13 )/2 c 13 /v^ 
-(l + s 13 )/2 c 13 /V2 J 



(Al) 



and 



( c 13 c u /V2 c 13 /y/2 si 3 

^ = -(1 + Si 3 )ci 4 /2 (l-Si 3 )/2 c 13 /V2 
(l-Si 3 )c 14 /2 _(l + Sl3 )/2 ci 3 /v^ 

\ -Sl4 

respectively, where we put the CP violating phase 5 i3 = [also 5 U = , as before in Eq. 
(2)]- 

In such a case, the neutrino oscillation formulae (9) lead for m\ ~ m| to the proba- 
bilities 



C13S14/V2 \ 

-(l + Si 3 )W 2 
(1 - s 13 )su/2 

Cl4 / 



(A 2) 



n/ \ 1 42-2 2 /ii2\/o 2-2 2 2 • 2 \ 

P(i/ e -> z/ e ) ~ 1 - c 13 c 14 sin x 2 i - c 13 (1 + c 14 J (2s 13 sin x 32 + c 13 s 14 sin x 4 i J 



P(v„ 



Ve) ^ 



O 2 2 2 -2 
— ZC 13 S 13 S 14 Sin X43 , 

1 1 4 2 • 2 

1 - - c 13 c 14 sin x 21 - 
1 



-(i + s y(i + 4) TSl3 4 



c 13 sin x 32 + -(l±s 13 ) s 14 sin x 4 i 



x 

1 4 2 • 2 2 

2 c i3 c i4 sm x 2 i + c 13 



-ic 2 3 (l±si 3 ) 2 s 2 4 sm 2 x 4 3 , 



«14 =F Sl3(l+C 2 



X 



=FSi 3 sin 2 x 32 + - ( 1 ± S13) s 2 4 sin 2 x 4 i 



-i^ 3 c 2 4 sin 2 X2i + i^ 3 (l r- 



± c 2 3 si 3 (l ± si 3 )s 2 4 sin 2 x 4 3 , 
1 



I'M sm 2 x 32 - -s 2 4 sin 2 x 41 



. 1 4 2 -2 
+ 2 C 13 S 14 Sm ^43 • 



(A3) 



The probabilities involving the sterile neutrino v s are 
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P(v»,t -> v a ) = (1 ± Si 3 ) 2 c 2 4 s 2 4 sin 2 x 4 i , 
P(z/ e ->• i/ s ) = 2c 2 3 c 2 4 s 2 4 sin 2 z 41 , 

P(i/ S -> = 1 - 4c 2 4 s 2 4 sin 2 x 4 i . (A 4) 

The formulae (A 3) and (A 4) are obviously reduced to Eqs. (11) and (12) in the limit of 
S13 -> 0. 

In the linear approximation with respect to a small S13 , only the probabilities P{y^ T — > 
^,t), P{ v h,t ~ ¥ v e) an d P{ v h,t ~ ¥ v s) of those given in Eqs. (A 3) and (A 4) get corrections 
to Eq. (11) and (12): 



5P(v^ T -> z/ ftT ) ~ ±si 3 Sm ( sin2 ^32 - c 2 4 sin 2 z 41 - sin 2 X43) , 
SP(u^ T -> i/ e ) ~ =FSi3S?4 ( sin2 ^32 + c 2 4 sin 2 x 4 i - sin 2 x 43 ) , 



5P{y^ T -> i/ s ) ~ ±2si 3 s 14 c 14 sin x 4 i . 



(A 5) 



Thus, in the linear approximation in S13, when making use of Eqs. (14), (29) and (18), 
(22) as well as (A 5), we obtain 



P(v e -> v e ) sol 
P(V P 



1 - c 14 sin (x 2 i)soi - ~(1 + c u) s i4 , 



^e)chooz ^ 1 - -(1 + C 2 4 )S 2 4 , 

P(z^ -> ^)atm ^ 1 - -(1 + c 2 4 ) sin 2 (a;32)atm - -(3 + c 2 4 )s 2 4 

z o 



+ s 13 s 14 
'I 



sin 2 (a;3 2 ) atm - -(1 + c 2 4 ) 



1 - 



1 + C 2 4 ) — Sl3-S 2 4 



Sin (X 3 2jatm 



^(3 + 4) + ^ 13 (l + c 2 4) 



s 2 

A 14 ) 



P(i/ jU -> z/ 6 )lsnd - -Si 4 sin 2 (a;4i)LSND - si3S 2 4 c i4 sin 2 (x 4 i)LSND - sin 2 (x 43 )LSND 



■4 2 2 

: s 14 S 13 S 14 C 14 



sin 2 (x 4 i) LSND + si 3 s 2 4 sin 2 (a; 4 3) L s ND , (a 6) 
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where Am^ < Am§ 2 < \Am 2 41 \ , |Am^ 3 | and (x 2 i) so i ~ 1, (^32)atm ~ 1, |(^4i)lsnd| ~ 1. 
We can see that, in this linear approximation, though the Chooz restriction (l + c 2 4 )s 2 4 ~ 
0.1 i.e., \s\ A ~ 1.3 x 10~ 3 , is not changed formally, the resulting amplitude for the possible 
LSND effect becomes slightly different. In fact, if (Am 2 ^] ~ |Am| 3 |, what holds both in 
the case of m\ <C m\ ~ m 2 , ~ m\ ~ 1 eV 2 and m\ ~ m\ <C m\ <C m\ ~ 1 eV 2 , the last 
Eq. (A 6) gives 

P(z/ M -> z/ e ) LS ND ^ ^ s i4(! + s i3) sin 2 (aj 4 i) lsnd , (A 7) 

leading to the oscillation amplitude 

sin 2 2# LSN D = ^ 4 (l + si 3 ) , (A 8) 

a little larger (smaller) than before with si 3 = 0, if S13 > (§13 < 0). 
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